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ABSTRACT 

LS 5039, a possible black hole x-ray binary, was recently observed with Giant Meter- 
wave Radio Telescope. The observed spectrum presented here shows that the spectrum 
is inverted at the low frequency. When combined with the archival data with orbital 
phase similar to the present observations, it shows a clear indication of a spectral 
turnover. The combined data are fitted with a broken power-law and the break fre- 
quency signifies a possible spectral turnover of the spectrum around 964 MHz. Truly 
simultaneous observations in radio wavelength covering a wide range of frequencies 
are required to fix the spectrum and the spectral turn over which will play a crucial 
role in developing a deeper understanding of the radio emitting jet in LS 5039. 
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1 INTRODUCTION 



LS5039 is a X-ray binary system consisting of a compact 
object and a companion star of ON6.5V(f) type with mass 
22.9M0. The mass of the compact object is not accurately 
known. But rece nt study shows that th e mass of the compact 
object is 3.7M0 (|Casares et al.l l|2005l ')'). Although t here are 
lot of debates regading the true nature of LS 5039 (jMirabeJ 
(2006 )). the dynamical mass estimate indicates that it is a 
, possible black hole x-ray binary can didate. The orbital p e- 
riod of the binary system is 3.9 days (jCasares et al.l ((20051)). 
The binary orbit is highly eccentric with an eccentricity of 
0.35 a nd has an inclination angle 24''.9±2''.8 (|Casares et al.l 

liooi)). 

LS5039 wa s proposed to b e a ve ry high energy gamma- 
ray emitter by IParedes et al.l (|200G| ) and was recently seen 
at TeV 7-ray energies by gamma-ray telesc o pe Hi gh En- 
ergy Stereoscopic System IjAharonian et al.l (120061 )). The 
gamma-ray light curve shows an orbital modulation of ~ 3.9 
days. The gamma-ray flux as well as the spectral hardness 
vary with the orbital phase. The flux variation with orbital 
phase can be accounted for by considering the absorption 
of gamma-rays due to 77 pair production. This also implies 
that 7-photons are emitted from a region within 1 AU of 
the compact object (jAharonian et al.l (|200y)). The observed 
spectrum is very steep at the superior conjunction and be- 
comes harder as the compact object moves away from it 
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and the spectral index goes below 2.0 as the compact object 
reaches the inferior conjuntion. But this variation of spectral 
hardness can not be adequately explained by the gamma-ray 
absorption via pair production alone. 

This source w as observ ed in X-ray by ROSAT 
Motch et al.l (I1997J)). R XTE (iRibo et all il99^) ASC A 



Mart occhia et al.' ('2OO5I )), Bepp o SAX pleis et al.1 l|2003l) ). 
Chandra ( Bosch-Ramon et al.l (|2005| )) and XMM 



Newton (|Martocchia et al.l ( 20051 )) missions during 1996 to 
2003. These observations were effectively carried out dur- 
ing different orbital phases of the binary system. All these 
observations revealed different flux levels and different spec- 
tral indices in the X-ray energy band, but no x-ray spec- 
tral state tra nsition were reported in th ese observations. 
More recently iBosch-Ramon et al.l (|2005l ) reported results 
of RXTE observations when the source was observed con- 
secutively for four days in July 2003. The X-ray flux is found 
to vary with the orbital phases and it maximizes near pe- 
riastron passage. This phase dependence of X-ray flux is 
due to the motion of the compact object accreting matter 
from the wind of the ma ssive companion wh il e mo ving in 
a highly eccentric orbit (JBosch-Ramon et al.l (|2005|)). But 
the observed anti-correlation between the photon index and 
the X-ray flux does not fit to the scenario where X-rays are 
considere d to be produced in and around an accretion disc. 
Bosch-Ra mon et al.l (120051 ) argued that the x-ray emission 
might be due to inverse Compton/synchrotron process in 
a relativistic jet which might possibly explain the observed 
anti-correlation between the photon index and the observed 
flux. But the photon indices reported in these observations 
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are very similar to that generally observed in case of hard 
state spectrum of black hole x-ra y binaries. 

LS5039 was first observed bv lMarti et~aLl (11998') in ra- 
dio with Very Large Array (VLA) and the observation re- 
sulted a power-law spectrum with a negative power-law in- 
dex. This indicates an optically thin synchrotron emission by 
non-thermal electrons. This observation a lso revealed a mod - 
erate variability in the radio flux. Later, iRibo et al.l (| 19991 ) 
carried out an observation campaign in radio with VLA at 
2.0, 3.6, 6.0 and 20.0 cm wavelengths and with Green Bank 
Interferometer (GBI) at 3.6 and 13.3 cm wavelengths. The 
observed spectrum was a power-law with a power-law in- 
dex of -0.46±0.01 supporting the previ ous observation of 
iMarti et al.] l|l998l 'l. IParedes et all (|2002h first resolved the 
source by observing with European VLBI Network (EVN) 
and the Multi-Element Radio-Linked Interferometer Net- 
work (MERLIN) and it was found that LS 5039 consists 
of an asymmetric two-sided jet w ith an maximum e xtension 
of one side ~ 1000 AU. Recentlv [Ribo et al.l (120081 ) studied 
the radio morphology of LS5039 using Very Long Baseline 
Array (VLBA) and it was found that the radio morphology 
changed into asymmetric from a very symmetric configura- 
tion within five days, but the radio emission from the core 
did not vary appreciably. But what is most interesting found 
m these radio observations is the optically thin non-thermal 
radio spectrum of the source which is in direct contradic- 
tion to the generally observed trend of flat or inverted radio 
spectr um for persiste n t black hole x - ray b inaries in the hard 
state \Fende\ 1(200 i ): \Gallo et al\ \200^) ). Therefore it is 
important to study LS5039 in the low-freguency radio band. 

We observed LS5039 at wavelengths 21 cm, 50 cm and 
128 cm (1280 MHz, 614 MHz and 234 MHz respectively) 
using Giant Meterwave Radio Telescope (GMRT). In this 
letter we discuss the results of the observation. In Section 
2 we discuss the observation and data reduction, results are 
discussed in Section 3 and we conclude the paper in Section 
4. 



2 OBSERVATIONS AND DATA REDUCTION 

The GMRT consist of 30 steerable antennas of 45m diameter 
in an approximate 'Y' shape similar to the VLA but with 
each antenna in a fixed position. Fourteen antennas are ran- 
domly placed within a central 1 km x 1 km square (the 'Cen- 
tral Square') and the remaining antennas form the irregular 
Y shape (six on each arm) over a total extent of about 25 
km. We refer the reader for detail s about the GMRT array 
to http://gmrt.ncra.tifr.res.in and lSwarup et al.l l|l99ll ). We 
observed LS 5039 at 1280 MHz on February 23, 2006 and 
simultaneously observed at 234 and 614 MHz on March 13, 
2006. 

The observations were made in standard fashion, with 
each source observation (30 min) interspersed with obser- 
vations of the phase calibrator (4 min) . The primary fiux 
density calibrator was either 3C 48 and 3C 286 , with fiux 
densities being on the scale of iBaars et al.l l|l977h . Either of 
these fiux calibrator was observed for 15 min at the begin- 
ning and end of each observing session. 

The data recorded with the GMRT was converted to 
FITS format and analyzed with the Astronomical Image 
Processing System (AIPS) using standard procedures. Dur- 



Table 1. Log of GMRT observations and results 



Frequency 

(MHz) 



Flux Density 
(mjy) 



MJD 



Phase 
range 



234 17.0±1.1 53788.78 0.61-0.67 

614 34.2±2.8 53788.78 0.61 - 0.67 

1280 38.1±2.3 53806.73 0.02 - 0.06 



ing the data reduction, we had to discard a fraction of the 
data affected by radio frequency interference and system 
malfunctions. After editing the data were calibrated and col- 
lapsed into a fewer channels. A self-calibration on the data 
was used to correct for the phase related errors and improve 
the image quality. 



3 RESULTS AND DISCUSSION 

LS5039 was observed simultaneously at 234 and 614 MHz 
during the orbital phase 0.61-0.67 whereas the observation 
at 1280 MHz was carried out during the orbital phase 0.02- 
0.06. The details of the observation are tabulated in Tabl e 
1. As previously r eported observations (jMarti et al.l l|l998l ): 
iRibo et al.l l|l999l ): |Paredes et all (|2002l )) show that the radio 
emission is core dominated, therefore the major contribution 
to the radio flux detected in present observations where the 
source could not be resolved, is mainly due to the emission 
from the core. Given that the source did not undergo any X- 
ray state transition over the period of observation concerned, 
it can be safely assumed that the variability of the source 
spanning this frequency range is not significant. The plot 
of the spectrum is shown in Figure 1. It is evident that the 
spectrum is inverted and it shows the trend of a possible turn 
over around 1000 MHz. To have a much clearer picture of the 
spectrum we plot the present data with the arch ival data. We 
have chosen the data reported bv iMarti et al.l (|l99S ) which 
was obtained during the orbital phase 0.63 similar to that 
of our second spell of observation. The combined spectrum 
is shown in Figure 2. The combined spectrum is fitted with 
a broken power-law given by 



F{u)=Fou~ 



for v < Vi 
for V > Ub 



where, Vb is the break frequency, ai and a2 are the power- 
law indices and _Fo is the nor malization. We hav e used 
Levenberg-Marquardt algorithm ([Press et al.l ([1993)) to fit 
the spectrum. The fitted parameters are Fo = 0.282 ±0.183, 
Qi = -0.749±0.111, a2 = 0.429±0.008 and i^b = 964±104 
MHz. This shows that the source has an inverted spectrum 
with a positive spectral index of 0.749; inverted spectrum 
extends upto the break frequency at 964 MHz and then the 
spectrum becomes a power-law with a spectral index —0.429. 
The power-law spectrum for frequencies greater than 964 
MHz indicates optically thin synchrotron spectrum by non- 
thermal electrons with spectral index 1.86. The break fre- 
quency at 964 MHz gives an indication that the spectrum 
has a possible turnover around that frequency. This needs 
to be confirmed by further truly simultaneous observations 
covering the complete frequency range. 
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V (MHz) 
Figure 1. Radio spectrum of LS 5039 as observed with GMRT 



The inverted spectrum indicates an optically thick self- 
absorbed synchrotron spectrum. In ideal condition the self- 
absorbed, optically thick synchrotron spectrum produced 
by relativistic electrons in a magnetized plasma is given 
by a power-law with index 5/2. But in case of a magne- 
tized plasma in a relativistic jet the physical condition of 
the plasma varies from point to point and it is very much 
dependent on the flow pattern and the shape of the jet. 
Therefore the resultant synchrotron spectrum from an un- 
resolved jet differs substantially from the ideal situation and 
the power-law index of the self-absorbed, optically thick syn- 
chrotron spectrum d iffers from 5/2. While the model by 
Blandford & Konigl (1979) or disk-jet symbiosis model by 
Falcke fc Biermann ( ^ 19951 explain the flat spectrum radio 



sources, iFalckd ( 19961 ) discussed an extended disk-jet sym- 
biosis model to explain the inverted spectra from Sgr A* 
and M81*. They extended the disk-jet symbiosis model by 
including a longitudinal pressure gradient which induces a 
moderate acceleration to a free conical jet along its flow di- 
rection. This longitudinal expansion leads to energy losses 
due to adiabatic expansion and makes the proper velocity 
of the jet dependent on its distance from the central engine. 
This model explain s the inverted spe c tra o f Sgr A* (JFalckd 
l| 19961 ')') and M81* iJFalcke fc Markofj ll200a')') and may also 
be applicable to radio cores in general (|Falckel (|l996r )). More- 
over, it predicts, for a maximal jet (where the internal power 
of the jet equals the kinetic power), the slope of the inverted 
spectra for different inclination angle of the source. As the 
microquasars are sca led down vers ions of the radio galaxies, 
the generic model bv lFalckd (|l996l ) can be used to model the 
inverted spectrum of LS5039 presented here. But a detailed 
theoretical modelling which is required to constrain the jet 
parameters and the radiative processes active in the source, 
is beyond the scope of the present paper. 



4 CONCLUSION 

Here, we report the low frequency radio spectrum of a X- 
ray binary LS5039 as observed with GMRT. The spectrum 
is inverted and shows an indication of spectral turnover. 




Figure 2. Broken power-law fitting of the combined data set 
(filled circle : present data, filled triangle : archival data from 
Marti et al. (1998)) 



With the assumption that the variability of radio emission 
from the source is not appreciable, present data are plot- 
ted with the archival data. The fitting of the data with a 
broken power-law reveals a possible indication of a spectral 
turnover at 964 MHz and an inverted spectrum with spec- 
tral index 0.749. This estimation of turnover frequency is 
indeed obtained using the archival VLA data IjMarti et al.l 
([1993)) with the assumption that the source is not variable. 
But this estimation may depend on the variability of the 
source. To have an idea about the sensitivity of the spectral 
break frequency on flux variability, we introduced a moder- 
ate (~ 10%) flux variation to VLA data and fitted the spec- 
trum. The spectral fitting shows that the break frequency 
i/b does not change appreciably and it is well within the 
uncertainities of the statistical fit. 

As LS 5039 is a possible black hole x-ray binary per- 
sistently in low/hard x-ray spectral state, the results of the 
present observation are consistent with the trend of inverted 
radio spectrum observed for persistent black hole x-ray bi- 
nary and this is the first source of this kind for which the 
indication of a spectral turnover is obtained. The spectral 
turnover can be used to determine the magnetic field in the 
jet provided the size of the source is known at that frequency. 
But it is to be noted that the magnetic field depends on 
the fourth power of the angular size. Therefore, a reliable 
measurement of magnetic field is possible if and only if the 
angular size of the source at the turnover frequency is deter- 
mined with sufficient accuracy. This result will further help 
to constrain the theoretical models to explain the broadband 
spectrum from the source. 

To have a deeper understanding on the radio spectrum 
and its dependence on the orbital phases it is important to 
have simultaneous observations over the entire radio band. 
Nevertheless, we are conducting a long term campaign to 
study the spectral behaviour at frequencies 234, 614 and 
1280 MHz over different orbital phases of LS5039. This wiU 
allow us to constrain the low frequency end of the spectrum 
of LS5039. 
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